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Abstract. A temporal analysis has been performed on a sample of 100 bright gamma-ray bursts (GRBs) with 
T9 o<2s from the BATSE Current Catalog. The GRBs were denoised using a median filter and subjected to an 
automated pulse selection algorithm as an objective way of identifing the effects of neighbouring pulses. The rise 
times, fall times, FWHM, pulse amplitudes and areas were measured and the frequency distributions are presented 
here. All are consistent with lognormal distributions. The distribution of time intervals between pulses is not 
random but consistent with a lognormal distribution. The time intervals between pulses and pulse amplitudes are 
highly correlated with each other. These results are in excellent agreement with a similar analysis that revealed 
lognormal distributions for pulse properties and correlated time intervals between pulses in bright GRBs with 
T9 o>2s. The two sub-classes of GRBs appear to have the same emission mechanism which is probably caused by 
internal shocks. They may not have the same progenitors because of the generic nature of the hreball model. 
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1. Introduction 


It has been recognised that GRBs may occur in two 
sub-classes based on spectral hardness and duration with 
T9 o>2s and T9o<2s (Kouveliotou et ah, 1993; Dezalay 


et ah, 1995; Paciesas et ah, 2001). The bimodal distribu¬ 
tion can be fit by two Gaussian distribu tions in the loga¬ 
rithmic durations ( McBreen et ah, 1994 ). There is signifi¬ 
cant evidence for a third subgroup as part of the long du¬ 
ration GRBs ( Mukherjee et ah, 1998 ; [Horvath, 1998 ) but 
this has been questioned because o f a possible BATSE se- 
l ection effect (Hakkila et ah, 2000 ). It has been suggested 
( jGline et ah, 1999 ) that the small group of GRBs with 
Tgo < 0.1 s form an additional ca tegory. The short GR Bs 
have a higher value of (U/14iax) ( Katz fc Canel, 1996 ), a 
much smaller value of the spectral lag ( Norris et ah, 2000|), 
a pulse shap e that depends on position in the burst ( jGuptq 
jet ah, 200C) and a smaller space density than long GRBs 
( ^chmidt, 2001 ). 

A variety of statistical methods have been applied to 
the temporal properties of GRBs with rgo>2s. It is im¬ 
portant to compare the temporal profiles of the long and 
short GRBs to determine the similarities and differences 
between the two classes in an objective way. A detailed 
objective analysis has been performed on the temporal 
profiles of a large sample of 319 bright GR Bs with T9o>2 s 
( iQuilligan et ah, 2001 ; Hurley et ah, 199^ ). The properties 
of the pulses in GRBs and the time intervals between them 
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were found to be consistent with lognormal distributions. 
These results can be used as templates for comparison 
with a similar analysis of GRBs with rgo<2s. 

The analysis method is presented in section 2 and the 
results in section 3. In section 4 the results are discussed 
and compared with the sample of long GRBs. 

2. Data Preparation 

The dataset was taken from the BATSE current catalog. 
The time tagged event data at 5 ms resolution was used. 
The data from the four energy channels were combined 
into a single channel to maximise the signal to noise ra¬ 
tio. A subset of the current catalogue was selected with 
GRB durations less than two seconds (Tgo < 2 s). These 
GRBs were ordered according to their peak flux in 256 ms 
(T 256 ms) and the first 100 bursts without data gaps and 
with 5 ms data available for the complete burst where se¬ 
lected . All GRBs in the sample had Tgsems >1.6 photons 

— 2 —1 
cm s . 

2.1. Background subtraction 

The first step in the data preparation involved selecting 
the appropriate background for subtraction. The start and 
end times for each burst were identified . A background 
section was then identified which was usually after the 
main section of the GRB. If 5 ms data was available for 
the background level estimation, it was used, otherwise a 
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Table 1. The parameters of the pulse properties in GRBs include the median value for the data, the median ^ and 
the standard deviation a expressed as natural logarithms, the width of the lognormal distributions at the ±50% level 
in normal space, the KS probability that properties of pulses selected at different levels of and Ti were drawn from 
the same distribution as those selected at > 5, r, > 50%. All pulses with > 5 were used for the time intervals 
and not just isolated pulses. The values for the pulse areas and amplitudes were obtained from 55 GRBs that were 
summed over two BATSE Large Area Detectors. 


Property 

Median 


cr 

Width (±50%) 

KS (r,, > 3,ri > 20%) 

KS (r, > 8, Ti > 80%) 

Rise Time 

0.035 

-3.31 

0.94 

0.012-0.11 

0.48 

0.16 

Fall Time 

0.056 

-2.89 

0.98 

0.017-0.176 

0.07 

0.19 

FWHM 

0.045 

-3.17 

1.01 

0.013-0.138 

0.73 

0.48 

Area 

1.5x10® 

11.7 

1.24 

28-520(xl0®) 

0.44 

0.15 

Pulse Amplitude 

1.02xl0'‘ 

9.22 

0.83 

3.8-27(xl0®) 

0.90 

0.72 

Asymmetry Ratio 

0.65 

-0.42 

0.91 

0.19-2.22 

0.78 

0.30 

AT 

0.095 

-2.24 

0.87 

0.038-0.3 

0.83 

0.96 


64 ms section was selected. The median value of the data 
points in the background section was subtracted from all 
other values in the active section of the GRB. 


2.2. Denoising and pulse selection 

A median hlter was used to denoise the GRBs. This de¬ 
noises a signal by finding the median of each bin. The best 
denoised signal was found by varying the bin size and the 
value of a. The wavelet method ( Quilligan et ah, 2001 ) 
was not used because the short durations of the GRBs 
produced a smaller number of data points than that re¬ 
quired for our wavelet method. 

The pulse selection method was described elsewhere 
( iQuilligan et al., 2001 ). The pulses selected had a thresh¬ 
old of 5 tJ (tct > 5) and isolated from adjacent pulses by at 
least 50% (ti > 50%). A value of ri > 50% implies that the 
two minima on either side of the pulse maximum must be 
at or below half the maximum value. A total of 313 pulses 
were selected with > 5 and 181 of these had r, > 50%. 


3. Results 

3.1. Distributions oft,-, tf, FWHM, pulse area, pulse 
amplitude and A T 

The distribution of the number of pulses with Tq- > 5 is 
given in Fig. 1. The median value of the number of pulses 
per GRB is 2.5 and is smaller than the value of 6 for 
long GRBs ( IQuilligan et al., 2001 ). The distributions of 


rise time t,, fall time tf and full width at half maximum 
(FWHM) for the isolated pulses are presented in Fig. 2 
along with the lognormal fits to the data. The distribu¬ 
tions of pulse amplitudes and areas are presented in Fig. 3 
for the isolated pulses. The distribution of time intervals 
between pulses with To- > 5 is given in Fig. 4. 

The median values and the parameters of the lognor¬ 
mal fits are presented in Table 1. The distribution of the 
pulse asymmetry ratios, t^/U , is not presented but the 
values are listed in Table 1. The pulse fall time is gener¬ 
ally longer than the rise time. The widths or ranges of the 
lognormal distributions at the ±50% are are also given in 


Table 1. The Kolmogorov-Smirnov (KS) test was used to 
indicate the effects of varying and t,. The pulse prop¬ 
erty distributions at To- > 5 and t, > 50% were compared 
with the distributions selected over large range of To- and 
Tj. The significance level probabilities from the KS test for 
the comparison are presented in Table 1 for To- > 3, t, > 
20 % and To- > 8, r, > 80% and all have acceptable values. 
It was found that when t, increased the pulses become less 
contaminated by overlaps from adjacent pulses. 


3.2. Correlation Analysis 

The correlation coefficients between a range of GRB pa¬ 
rameters and pulse parameters were obtained using the 
Spearman Rank Order correlation coefhcents p with as¬ 
sociated probabilities and the values are listed in Table 
2. The correlation coefficients for the time intervals be¬ 
tween pulses are listed in Table 3 and they are strongly 
correlated with each other. 
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Fig. 1. The number of pulses versus number of y-ray 
bursts. 
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Fig. 2. The distribution of U, <f, and FWHM (a-c) for 
pulses with > 5 and n > 50%. The dashed lines are 
lognormals that are consistent with the data. 


4. Discussion 


The major result is that pulses in GRBs with Tgo<2 s have 
very similar distributions of tr, U, FWHM, pulse ampli¬ 
tude, pulse area and AT. The frequency distributions are 
broad and all are compatible with lognormal distributions, 
and at the 50% level span an order of magnitude in the 
particular pulse property (Table 1). The same pulse prop¬ 
erties in GRBs with T9 o>2 s are also well des cribed by log¬ 
normal distributions ( Quilligan et ah, 2001 ). The results 
of the correlation analysis on the GRBs and the properties 
of the pulses (Table 2) are generally in good agreement be¬ 
tween long and short GRBS. The main difference is the 


Table 2. Spearman Rank Order correlation coefficients 
p and associated probabilities between a range of burst 
(first 3 rows) and pulse parameters (second 10 rows). 


Properties 

P 

Probability 

No. of Pulses vs. Tgo 

0.18 

0.073 

No. of Pulses vs. Total Fluence 

0.38 

9.8 X 10"® 

Tgo vs Total Fluence 

0.21 

0.037 

Rise Time vs. Fall Time 

0.54 

2.7 X lO”^'" 

Rise Time vs. FWHM 

0.79 

1.3 X 10-'‘° 

Rise Time vs. Area 

0.59 

3.3 X 10"“ 

Rise Time vs. Pulse Amplitude 

0.013 

0.90 

Fall Time vs. FWHM 

0.75 

2.6 X 10"®® 

Fall Time vs. Area 

0.64 

2.0 X 10"^^ 

Fall Time vs. Pulse Amplitude 

0.078 

0.45 

FWHM vs. Area 

0.60 

6.89 X 10"“ 

FWHM vs. Pulse Amplitude 

-0.05 

0.64 

Area vs. Pulse Amplitude 

0.65 

1.1 X 10"^^ 


Table 3. The Spearman correlation coefficients p and 
probabilities for ATs (first 2 rows) and PAs (second 2 
rows). The first and third rows refer to adjacent ATs and 
PAs while the second and fourth rows refer to ATs or PAs 
separated by one time interval or pulse. The first / second 
entry for the ATs gives p when the ATs are unnormalised 
/ normalised by Tgo. The first / second entry for the PAs 
gives p when the PAs are unnormalised / normalised by 
the maximum peak in the burst. 


Total Number 

P 

Probability 

140 

0.42/0.30 

1.7 X 10"^V3-0 X 10"® 

84 

0.48/0.32 

4.9 X 10"®/3.2 X 10"® 

213/107 

0.39/0.24 

5 X 10"®/1.1 X 10"^ 

140/84 

0.13/-0.09 

0.13/0.42 


pulse amplitude because it is not anticorrelated with tr, 
tf and FWHM. In GRBs with Tgo>2s the values of AT 
are not random but consistent with a lognormal distribu¬ 
tion with a Pareto-Levy tail for a small number of long 
time intervals in excess of 15 s. It was also found that the 
values of AT were correlated over most of the GRBs and 
up to more than 20 pulses in the GRBs with large num¬ 
bers of pulses. In short GRBs there is a strong correlation 
between the values of AT and also between the pulse am¬ 
plitudes (Table 3) but the latter are better correlated in 
GRBs with T9o>2 s. There are remarkable similarities be¬ 
tween the statistical properties of the two classes of GRBs. 
The clear conclusion is that the same emission mechanism 
accounts for the two types of GRBs. This conclusion is in 
agreement with a very different analysis of the temporal 
structure of short GRBs (Nakar & Piran, 2001). The ex¬ 
ternal shock model (Dermer & Mitman, 1999) has serious 
difficulties in accounting for GRBs with Tgo <2 s and with 
the non-random distribution of correlated time intervals 
between pulses. The results presented here provi de con- 
siderable support for the internal shock model (Rees fc 
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et al., 2001). The similar distributions of AT in short and 
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Fig. 3. The distributions of pulse area and amplitude 
selected at > 5 and ri > 50% and consistent lognormals 
(dashed lines). These distributions were generated from 
55 GRBs that were summed over two BATSE Large Area 
Detectors. 


long GRBs are indirect evidence for rotation powered sys¬ 
tems with super-strong magnectic fields. 

There is strong evidence that long GRBs originate in 
star forming galaxies providing evidence that they are 
linked to massive stars and supernovae (Meszaros, 2001). 
A variety of models, such as collapsars and hypernovae 
have been proposed as progenitors (e.g. MacFadyen & 
Woosley 1999). Models of these systems imply that jets 
can be generated to produce GRBs by an internal shock 
model. These models cannot account for GRBs with Tgo 
less than about 5 s. The shorter GRBs have been at¬ 
tributed to neutron star (NS) - NS mergers or NS - black 
hole mergers (e.g. Ruffert & Janka 1999). A common in¬ 
gredient of most GRB models is the formation of a black 
hole with a temporary torus that accretes and powers the 
relativistic jets and the GRBs. The overall energetics of 
the various progenitor models differ by at most an order of 
magnitude. The GRBs appear to have the same emission 
mechanism and possibly different progenitors for long and 
short bursts. 

5. Conclusions 

A sample of bright GRBs with Tgo <2 s have been denoised 
and analysed by an automatic pulse selection algorithm. 
The results show that the distribution of the properties 
of isolated pulses and time intervals between all pulses 
are compatible with lognormal distributions. The same 
mechanism seems to be responsible for both long and short 
GRBs and is attributed to the internal shock model. 
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